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ABSTRACT 

We study the temperature structure of the rich cluster of galaxies Abell 1367 using two independent 
methods of correcting for the ASCA PSF. The results of the two methods are in excellent agreement and 
give solid evidence of a strong and localized shock in the intracluster medium (ICM). Our analysis sug- 
gests that we are observing the merger of two subclusters. We find that the lower luminosity subcluster, 
located to the northwest, has a higher temperature of 4.2 ± 0.3 keV, compared to the more luminous 
southeast subcluster, whose temperature is 3.2 ± 0.1 keV. An analysis of the ROSAT surface brightness 
profiles of both subclusters is also presented. The data agree well with predictions from numerical simu- 
lations of mergers between subcluster sized masses, in particular, with early stages of the merger, prior 
to the first core passage. 

Subject headings: galaxies: clusters: individual (A1367) — intergalactic medium — X-rays: galaxies — 
methods: data analysis 



1. INTRODUCTION 

Optical and X-ray studies (e.g. Geller & Beers 1982; 
Forman et al. 1981; Dressier & Shectman 1988; Jones 
& Forman 1992; Mohr et al. 1994; Bird 1994; Slezak et 
al. 1994) have shown that galaxy clusters are dynami- 
cally evolving systems exhibiting a variety of substructure. 
Thus, we expect to see key indentifiers of the merging pro- 
cess in the temperature maps of galaxy clusters with sig- 
nificant substructure. In particular, in the early stages 
of an unequal merger, i.e. one subcluster larger than the 
other, simulations show significant heating of the smaller 
subcluster above the temperature of the local ICM (Evrard 
1990a and b; Schindler & Miiller 1993) as well as the de- 
velopment of a shock located between the two subclusters. 

Abell 1367 has previously been identified as having sig- 
nificant optical and X-ray substructure (Bechtold et al. 
1983; Grebenev et al. 1995). The X-ray emission is 
elongated along a southeast-northwest axis, and contains 
small, localized "clumps". The cluster has a relatively 
cool gas temperature and a high spiral fraction (see Bah- 
call 1977 and Forman & Jones 1982) typical of what is 
expected for a dynamically young system. 

In this paper we report on the analysis of the struc- 
ture of the A1367 galaxy cluster as mapped by the X-ray 
emission observed with the ROSAT PSPC and the ASCA 
detectors. In Section y, we apply two independent meth- 
ods which account for the energy dependent ASCA Point- 
Spread-Function(PSF) and produce moderate spatial res- 
olution {^ 4') temperature maps. This is an extension 
of the work presented by Churazov et al. (1996a). We 
also describe our surface brightness analysis of the ROSAT 
PSPC data. Section g presents the results of the temper- 
ature determinations and our estimates of the masses of 
each subcluster. Finally, Section || briefly summarizes our 



results. All distance dependent quantities have assumed 
Ho — 50 km s~'^ Mpc~^, qo = 0.5, all coordinates are given 
in the J2000 system, and unless otherwise noted all error 
bars are 67% confidence level (Icr) errors. 

2. OBSERVATIONS & METHODS 
2.1. ROSAT Analysis 

The ROSAT PSPC observed A1367 from 29 Novem- 
ber to 2 December 1991 (RP800153); observing details are 
listed in Table |l|. We corrected the PSPC image for tele- 
scope vignetting as well as the removal of times with high 
solar/particle backgrounds using the standard procedures 
outlined by Snowden (1994; also Snowden et al. 1994). By 
combining only the data from Snowden bands 4 through 
7 (0.44-2.04 keV), we excluded the lower energies, which 
are likely to have higher background X-ray contamination. 
The background was not subtracted during this stage of 
the processing; instead, it was included as a constant com- 
ponent in our fits of the surface brightness, as described 
below. The central area of the PSPC image is shown in 
Figure 0, after having been smoothed with a 30" Gaussian. 

The extended cluster emission appears peaked in two lo- 
cations, a primary peak near ll'' 44.8™ -|-19°42'(hereafter 
the SE subcluster) and a secondary peak towards the 
northwest at ll'' 44.4™ -Hl9°52'(hereafter the NW sub- 
cluster). The temperature distribution and previous radio 
results- both discussed later in this paper- suggest that 
we are observing the merger of two subclusters. Because 
of this, we performed a surface brightness analysis deriv- 
ing the X-ray surface brightness profile and the core radius 
assuming a standard /3-model on each region independent 
of and excluding the other. This data also were used as a 
surface brightness model for one of the spectroscopic anal- 
yses of the ASCA data (Method B). 
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To aid in the location of potential point sources contam- 
inating the field, we smoothed the PSPC image on scales 
from 30" to 8'. The smaller scales were used to locate point 
sources near the center of the image, while the larger scales 
were used near the edges where the distortions of the PSF 
are large. Potential point sources were identified by eye 
and these areas were excluded from further analysis of the 
surface brightness distribution. 

To accurately locate the peaks of the two subclusters, 
we masked off one peak as well as the point sources, and 
centroided the other. We then repeated the process for 
the other peak. The equatorial coordinates (J2000) of the 
two peaks in the ROSAT image coordinate frame are: SE- 
a = ll''44"50" 5 = 19°41'44", NW- a = ll''44"22'* 6 = 
19°52'27". These appear to be offset +0.175'* in RA and 
+ 11" in Dec from the true sky position based upon the lo- 
cation of NGC 3862 (the bright X-ray point source in the 
southeast). This is roughly consistent with typical point- 
ing errors for ROSAT. 

To generate radial profiles, we defined 1' wide annuli, 
with inner radii from 0' to 46', centered on each peak. 
To avoid contamination of one subcluster by the other, 
we excluded the third of each annulus which was on the 
side towards the other subcluster. For the SE subcluster 
the excluded azimuths ranged from 255°to 15°and for the 
NW subcluster from 75°to 195°(with the angle measured 
counter-clockwise from North). 

We then measured the average surface brightness 
(cts s~^ arcmin^^) in each annulus, and fit the resultant 
surface brightness profiles with a standard hydrostatic, 
isothermal /3-model: 



E(r) = S„ 



' + '«: 



n-(3/3-^) 



(1) 



(Cavaliere & Fusco-Femiano 1976). Because of Al367's 
low redshift [z = 0.0215) the cluster emission overfills 
even the PSPC's large field of view. This makes direct 
measurement of the background difficult. Therefore, we 
included a constant background component in our model. 
The best fit backgrounds were (1.79 ± 0.17) x lO"'^ and 
(1.29 ± 0.29) X 10-* cts s-i arcmin'^ for the SE and 
NW subclusters respectively, which are both consistent 
with typical PSPC backgrounds. While the best fit back- 
ground for the SE subcluster is slightly larger, possibly 
indicating a small amount of contamination by the NW 
subcluster, the errors are consistent with a single constant 
background. The results of our fitting are given in Table 
and the fits themselves are shown in Figure 0. 

2.2. ASCA Analysis 

To correctly characterize the temperature from ASCA 
data at some location for an extended source, we must ac- 
count for the extended and energy dependent PSF of the 
telescope (Takahashi et al. 1995). The primary difficulty 
for spatially resolved spectroscopy is caused by the outer 
part, or wings, of the PSF. Failing to correct for this com- 
ponent can lead to spurious temperature and abundance 
gradients, although for A1367, due to its relatively low 
temperature and large spatial extext, these efi'ects are not 
expected to be very strong. 

We have employed two independent methods which 



account for the broad energy dependent ASCA Point- 
Spread-Function (PSF) to construct temperature maps for 
A1367. Method A (Churazov et al. 1997; Gilfanov et al. 
1997), provides a rapid approximate correction of the ex- 
tended wings of the PSF. In contrast. Method B (Marke- 
vitch et al. 1997) performs an exact convolution of a sur- 
face brightness model with the PSF and the effective area 
of the telescope to generate model spectra which are then 
compared to the data. 

ASCA observed A1367 on 4-5 December 1993, with four 
pointings, one pair centered on the northwest region and 
the other centered on the southeast. Within each pair, the 
two pointings were offset from each other by 1.89' along a 
roughly SE-NW axis. This offset allowed an evaluation of 
any systematic effects. Details of the ASCA observations 
also are given in Table |l|. 

Preliminary to correcting the PSF, the ASCA data were 
"cleaned" with standard processing tools (Arnaud 1993). 
A cutoff rigidity of 8 GeV/c, minimum Earth elevation an- 
gles of 5° for the CIS and 20° for the SIS, and maximum 
count rate of 50 cts/s in the radiation belt monitor were 
used. The GIS's background was generated from an ap- 
propriately weighted combination of background maps for 
all rigidities, and for the SIS's, all hot-pixel events were 
removed and the background maps were scaled by total 
exposure. 

2.2.1. Approximate Fitting of the Wings of the ASCA 
PSF- Method A 

Method A approximates the ASCA PSF as having a 
core and broader wings (see Churazov et al. 1997 and 
Gilfanov et al. 1997 for details). The core PSF is cor- 
rected explicitly while the energy dependant PSF of the 
wings is corrected using a Monte Carlo algorithm. Af- 
ter the approximate subtraction of the scattered flux in 
the wings of the PSF, the temperature is determined us- 
ing one of two approaches. The first fits the spectrum in 
each 15" pixel with a linear combination of two fiducial 
single-temperature spectra and then smooths the result to 
reduce the noise (Churazov et al. 1996b). This approach 
yields a continuous (unbinned) temperature map of the 
cluster. Central to this method is the fact that thermal 
spectra having typical cluster temperatures (^ 2 keV) can 
be approximated as a linear combination of two spectra 
bounding the temperature range in the cluster (Churazov 
et al. 1997). For A1367 we used fiducial spectra with 
kT=2 and 6 keV. The results for the CIS data for A1367 
are shown in Figure 0, along with intensity contours from 
the ASCA data. Only temperatures with a ^^^^ < 1.5 
are shown. Results for the SIS are similar, although with 
a smaller detector field of view. 

The second temperature fitting approach proceeds by 
defining a series of regions which we chose to lie strad- 
dling the line connecting the two subclusters (see Fig- 
ure 0). Each region was 5' in width and approximately 16' 
in length. We excluded a region 5' in diameter around the 
bright point source (NGC 3862) near the SE subcluster, 
so as to prevent its signal from contaminating our temper- 
ature fits. This diameter was chosen to be large enough 
to more than contain the central part of the PSF. In ad- 
dition, we estimate that the point source contributes only 
about 33% of the total flux in this excluded region, which 
gives us confidence that there is little contamination in the 
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neighboring regions (#1 and #2). All of the data in each 
region were binned together, and then output as compos- 
ite spectra. These spectra were then fit with a standard 
Raymond-Smith model using the XSPEC package utiliz- 
ing the data from 1.5-2.0 and 2.5-11.0 keV. These limits 
were chosen to exclude the poorly calibrated region near 
the gold edge at 2.2 keV, and the extreme low energies for 
which the PSF remains poorly determined. 

2.2.2. Multiple Region Simultaneous Spectral Fitting- 
Method B 

Method B employs a model surface brightness distribu- 
tion which is convolved with the mirror effective area and 
the ASCA PSF (Takahashi et al. 1995) to produce model 
spectra for a set of user defined regions in the ASCA detec- 
tor planes. The spectra from the desired regions arc then 
fit simultaneously. A more detailed description, including 
a discussion of systematic uncertainties, can be found in 
Markevitch et al. (1996 and 1997). 

Integral to this method is the use of a detailed surface 
brightness model. To this end we used the ROSAT PSPC 
image from our surface brightness analysis, blocked, ro- 
tated and shifted to coincide with the ASCA image. We 
chose the regions used in this analysis to be identical to 
those defined in Method A. 

To perform the actual temperature fitting, the pulse 
height data were binned in energy to achieve an adequate 
signal to noise ratio. The same energy range used in 
Method A was utilized, with the bins defined to be 1.5- 
2.0, 2.5-3.5, 3.5-5.5 and 5.5-11.0 keV. AU of the process- 
ing so far was performed independently for each detector 
(SIS-0, SIS-1, GIS-2, and GIS-3) for each pointing. To de- 
termine the temperatures in each of the model regions, all 
of the SIS data- from both the and 1 detectors- for all 
four pointings, for all the regions and for all of the energy 
bands were simultaneously fit . A similar procedure was 
followed for the CIS. 

To estimate the errors in our temperature solutions, we 
measured the standard deviation of the distribution of fit 
temperatures from 200 simulated spectra. Each spectrum 
was constructed by performing a Monte-Carlo simulation 
of the counts in each energy band, assuming a Gaussian 
distribution about the observed number of counts in that 
energy band in the data. The systematic errors were added 
to the data and model as appropriate. 

3. RESULTS & DISCUSSION 
3.1. Temperature Structure 

Figure H shows the continuous temperature map pro- 
duced by Method A and indicates a temperature gradient 
across the cluster, increasing from the 3.0 keV in the SE 
to 4.3 keV in the NW. The regions defined earlier were de- 
signed to assess the significance of this trend, spanning the 
cluster along an axis from the SE to the NW (see Figure |3 
and Figure ph . Both spectral analysis methods (A & B) 
discussed above were applied to the data and the results 
of the fitting are given in Table || and Figure |. 

The first feature of note in Figure ^ is the excellent 
agreement between the two methods for both the GIS and 
SIS detectors. This gives us confidence that the results of 
our fitting procedure are correct, at least to our knowledge 
of the ASCA PSF. 



The next feature to note in Figure ^ is that the temper- 
ature variation appears to be abrupt rather than smooth, 
with a jump occuring in region #4. The spectral fits for 
the SE subcluster (regions #1-3) are consistent with a con- 
stant temperature of 3.2 ±0.1 keV, and although the data 
for the NW subcluster (regions #5-7) have larger error 
bars, they are consistent with a constant, but higher, tem- 
perature of 4.2±0.3 keV, spanning the bulk of the emission. 

To study the nature of the transition and as a consis- 
tency check, we shifted the regions by one half of a box 
width (2.5') along the SE-NW axis and reapphed Method 
A to the GIS data. For the shifted regions that overlapped 
with the original regions #l-#3, i.e. the SE subcluster, 
the fit temperatures were unchanged. Similarly the tem- 
peratures for the NW subcluster, regions #5-#7, were un- 
changed. However, the abruptness of the transition from 
SE to NW became more pronounced. 

Focusing on the regions near the transition- #3, #4 and 
#5- the original GIS Method A temperatures were 3.3, 3.6 
and 3.8 keV respectively. After the shift, the region to the 
southeast of the middle of #4 had a temperature of 3.4 
keV while the temperature of the region to the northwest 
was 3.9 keV. 

If the extent of the transition was large, the shifted re- 
gions would, by their partial inclusion of intermediate gas 
in two shifted regions, have had a higher temperature than 
before on the southeastern side and a lower temperature 
on the northwestern side. Instead the intermediate region 
effectively disappeared. Region #4 appears to have been 
intermediate because it contained nearly equal amounts of 
the SE and NW subclusters. 

This abrupt temperature change indicated by the tem- 
perature fits with and without the shift is strongly sug- 
gestive of a shock located nearly at the midpoint of region 
#4, that has been generated during a collision between the 
two galaxy subclusters. In fact the temperature distribu- 
tion and intensity contours are very similar in nature to 
cross-sectional temperature and projected surface bright- 
ness profiles shown in Figures 3b-c and 5b-c respectively 
in Schindler & Miiller (1993). In these simulations, at 0.95 
Gyr after the beginning of the merger, the smaller subclus- 
ter shows extensive heating in its core relative to its initial 
state as well as a strong gradient with radius. By 2.66 Gyr 
the cores are nearly in contact and the thermal gradient in 
the smaller subcluster has dissipated considerably. At the 
same time the temperature of the core of larger member 
of the merger has begun to increase and develop a thermal 
gradient. 

In A1367, the measured gas temperatures of both the SE 
and NW subclusters are hotter than that expected from 
the luminosity-temperature relation for clusters (e.g. Edge 
& Stewart 1991; David et al. 1993). Although there is 
moderate dispersion in the Lx-T relation it is possible that 
the SE subcluster gas also has begun to be heated due to 
the merger. This potential heating of the SE subcluster 
and apparent lack of thermal gradients in both subclusters 
as well as the clear separation of the two peaks suggests 
that we are observing the merger at a stage intermediate 
to those shown in the figures of Schindler & Miiller, at 
approximately 1.8 Gyr after the onset. 

Observations by Gavazzi et al. (1995) also suggest that 
A1367 is currently undergoing a merger. They find three 
head-tail radio galaxies in the NW subcluster that have 
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extremely large radio/IR flux ratios as well as an extreme 
excess of giant HIT regions on their leading edges, all of 
which are pointed towards the SE subcluster. From sim- 
ulations, Roettiger et al. (1996) find that the expected 
shock generated from a merger event could "induce a burst 
of star formation" as well as help to generate head-tail 
morphology. 

Finally, in an effort to generate a detailed map of the 
cluster we divided each of the rectangular regions into 
three roughly square (5'x5'.3) subregions, and re-applied 
Method A. By combining the temperature fitting for the 
GIS and SIS data we attempted to improve the statistics in 
each region and produce composite temperatures. The re- 
sults are shown in Figure ^ We also added four regions to 
the northeast of regions #l-#4 to examine the cool feature 
that appears in Figure g. We see a trend of decreasing tem- 
perature from west to east across the subregions around 
the SE subcluster (#41 to #14). This may indicate that 
the merger is slightly oblique rather than head-on; how- 
ever, we note that the variations are not highly significant. 

Figure || also helps to rule out the possibility that we 
are viewing the NW subcluster through a hot isothermal 
shell. In that case, we would expect the temperature to 
decline with decreasing projected radius, which is opposite 
to the trend found in our temperature maps. 

3.2. Density Profiles 

Typical values of the core radius for relaxed clusters 
range up to 0.6 Mpc with the peak of the distribution 
around 0.2 Mpc (Jones & Forman 1984). The core radii 
of the SE and NW subclusters are 0.42 and 0.49 Mpc re- 
spectively, and lie significantly toward the high end of the 
distribution. Simulations by Roettiger et al. (1996) show 
that during a merger the core radius of the gas can in- 
crease by a factor of two or more due to the increase in 
the central entropy of the gas through shocks. Further, 
the core radius of the NW subcluster is larger than that 
of the SE subcluster indicating that the NW subcluster is 
even farther from a relaxed state than the SE subcluster. 

Similarly, the best fit values of /3 are 0.73 and 0.66 for 
the SE and NW respectively, which are relatively large 
compared to relaxed clusters with temperatures similar to 
those of A1367. In a relaxed cluster /3 generally lies in the 
range of 0.4 to 0.8 and increases with the gas temperature. 
The typical value for a 4 keV cluster is between 0.5 and 
0.6 (Jones & Forman 1997). 

3.3. Abundances 

For Method A, where we fit the spectrum for each re- 
gion independently, we allowed the abundance to be a free 
parameter in the fit. The results with error bars are given 
in Table and in Figure o. Even though the error bars 
are quite large, there is some suggestion from these re- 
sults that the NW subcluster has a lower abundance than 
the SE subcluster. To test this we defined two regions on 
the GIS data, one surrounding each subcluster, and again 
applied Method A. The results from these fits gave abun- 
dances relative to solar of 0.26 ± 0.06 and 0.11 ± 0.05 for 
the SE and NW subclusters respectively. We note, how- 
ever, that large changes in the abundance, e.g. assuming 
an abundance of 0.3 solar for the NW subcluster, had no 
significant effect on the fit temperatures. 



In Method B, where the abundances were held fixed dur- 
ing the fits, we determined the temperatures twice, once 
with an abundance of 0.3 solar and once with an abun- 
dance of 0.4 solar. The differences in the temperature 
solutions were inconsequential. Because the temperatures 
are relatively insensitive to abundance, for Method B we 
have presented only the results using the typical abun- 
dance of 0.3 solar. 

3.4. Mass Estimates 

We made two estimates for the mass of each subcluster 
by applying the equation of hydrostatic equilibrium with 
the same density profiles, but two independent tempera- 
tures. In the case of a merger and the subsequent heating 
of the ICM gas, the mass estimates though very uncertain 
do provide a means of comparing the masses of the two 
subclusters. The first estimate uses the temperature we 
derived from our spectral fitting; the second came from 
applying the observed relation between temperature and 
X-ray luminosity found by David et al. (1993). Indepen- 
dent of these calculations we also estimated the masses of 
the X-ray emitting gas in each subcluster. 

The assumption of hydrostatic equilibrium and spheri- 
cal symmetry gives a simple equation for the total mass of 
the emitting system which when combined with a density 
profile given by an isothermal /3-model reduces to: 

Mir) = 1.13 X 1015/3 (^—] (^—] ^TZM^Mq. 
^' ^ \10keVj \MpcJ l + {r/Rc)^ " 

(2) 
To test the assertion of isothermality, we constructed 
three semi-circular annuli- radii 0-6', 6-12', and 12-18'- 
centered on each subcluster such that each sampled the 
region away from the other subcluster in order to exclude 
the region of the shock and the other subcluster- and ap- 
plied Method A to the data. Due to the limited coverage 
of the SIS in the outer areas, we used only the GIS data. 
A region 5' in diameter centered on NGC 3862 was once 
again excluded from the analysis. 

The summed spectra for each region were fit with a 
Raymond-Smith model to determine the temperatures. 
For the SE subcluster annular regions we find tempera- 
tures from Method A for the GIS data only of 3.1 ± 0.1 
keV, 3.1 ±0.2 keV and 3.6±0.4keV for the inner, interme- 
diate and outer annuli respectively. For the NW subclus- 
ter the Method A GIS only temperatures were 4.9 ± 0.5 
keV, 4.7 ± 0.4 keV and 4.6 ± 0.7 keV for the inner, in- 
termediate and outer annuli respectively. We note that 
the large errors in the outermost annuli are due to their 
close proximity to the edge of the GIS imaging area. Both 
profiles are consistent with isothermal conditions and with 
the previous measured temperatures from the appropriate 
rectangular regions. 

We know from our previous results that there is temper- 
ature structure in the cluster, namely the shock. However, 
to derive an estimate of the masses, we assume isother- 
mality throughout each subcluster. We used the weighted 
average temperatures (Tc) of 3.2 keV and 4.2 keV applied 
to Equation | to calculate the mass within 0.5 Mpc. These 
include the results from both Method A and Method B for 
the GIS and SIS detectors. We also have extrapolated the 
observed mass to that within 1 Mpc for ease of comparison 
with previous work. The temperature of 4.2 keV-instead 
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of those listed above- was selected for the NW subcluster 
due to the consistently higher results found for it (regions 
#6, #7 and #8) using Method A with the GIS data (see 
Figure U) . The resultant mass estimates are given in Ta- 
ble | 

The second hydrostratic equilibrium mass estimate uses 
the same gas density profile as the first, but with a different 
temperature. We have used the empirical relation between 
temperature and luminosity from David et al. (1993), 



kT, 



eff 



10 



-0.72 



/ Lbol \ 



0.297 



(3) 



to estimate an effective temperature (T^ff) for a similarly 
luminous, but undisturbed, cluster. 

We determined the total flux within the ROSAT band- 
from 0.5 to 2.0 keV- from each of the previously defined 
annuli. Applying azimuthal symmetry, we found the total 
flux within a radius of 18' (0.66 Mpc) for each subcluster, 
and then calculated the bolometric luminosity (David et 
al. 1997). 

Finally, the gas mass was estimated by inverting the 
formula from David et al. (1990), 



27rneni/Aoa3 
L(r) = TZ :r-r. X 



(1 - 3/3) 



-3/3+1 



(1 



-3/3+1 



ds (4) 



and solving for the central density. We then integrated the 
density distribution to find the gas mass. 



3/3 



Mgas = 47rpo 



Re 



ds, 



(5) 



using the density distribution corresponding to the surface 
brightness given in Equation n^ for an isothermal gas with 
Po = P'eTT'emp wherc fie is the mean molecular weight per 
electron. 

The resultant gas masses are given in Table and yield, 
when compared to the mass estimates (using Tcff for the 
NW subcluster), gas mass fractions of ^12% at 0.5 Mpc 
and ~16% at 1 Mpc in each subcluster which is typical of 
rich clusters. 

Although, as stated above, there is considerable disper- 
sion in the Lx — T relation, temperature estimates from the 
luminosity for both subclusters are lower than the mea- 
sured values. For the SE subcluster the difference is mod- 
erate with an estimated temperature- and thus mass- 31% 
lower than the actual measured values. However, for the 



NW subcluster the difference is nearly twice as large. We 
measure a luminosity of 0.29 x lO^'* ergs s~^ and thus esti- 
mate a gas temperature of 2.0 keV which is 52% lower than 
the measured value of 4.2 keV. Again the results for the 
masses within 0.5 Mpc and 1.0 Mpc are given in Table ||. 
As discussed above the estimates of the masses of the 
two subclusters have considerable uncertainty associated 
with them due to the perturbations from the merger. How- 
ever, the two subclusters appear to be experiencing similar 
changes to their density profiles, thus allowing a relative 
comparison of their masses. While the masses derived from 
the measured temperatures are essentially equal, the NW 
subcluster appears to have been heated more than the SE 
subcluster, thus suggesting that it is actually less massive 
than the SE subcluster. This is further supported by the 
calculated gas mass. 

4. CONCLUSION 

We have analyzed the ROSAT PSPC and ASCA SIS and 
GIS observations of A1367. For the ASCA data, we have 
applied two different analysis techniques for measuring the 
intracluster gas temperature and find excellent agreement. 

Our analysis indicates that we are observing the early 
stages of a slightly unequal merger between the two sub- 
clusters occuring along a SE-NW axis nearly in the plane of 
the sky. We find evidence for a shock-like feature along the 
merger axis between the two subclusters as well as heating 
of the gas throughout both subclusters, with the smaller 
NW subcluster being heated more than the SE subcluster. 
This is in excellent agreement with predictions from the 
merger simulations by Evrard (1990a and b) and Schindler 
& Miiller (1993). We also note that the surface brightness 
profile of the NW subcluster is "puffing out" as indicated 
by its larger core radius. This is similar to the effects iden- 
tified by Roettiger et al. (1996) as a signature of a merger 
event. 

Our detailed temperature map of the cluster suggests 
that the merger may be occuring slightly obliquely, with 
the cluster cores passing each other traveling north-south, 
but the statistics of the data are not sufficient to support 
any firm conclusion on this point. 

Future studies of A1367 and other merging clusters will 
provide a clearer picture of the detailed interactions which 
occur as clusters form. With the launch of AXAF in 1998, 
we can expect to obtain much higher angular resolution 
temperature maps to study the merger process and the 
structure of the shocks which are produced. 
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from the Smithsonian Institute and NASA contract NAS8- 
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Fig. 1. — Central l°xl°of the PSPC image of A1367 with no background subtraction. This scale was chosen to allow easy comparison 
with the ASCA data. The outline of the regions used in the spectral fitting for both methods are shown. Also shown is the excluded circular 
region surrounding NGC 3862. 
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Fig. 2. — /3-modcl fits (dashed fines) with eonstant baekgrounds (dash-dotted fines) for the two subclusters. Core radii are given in 
arcrninutes and are 0.49 and 0.42 Mpc for the NW and SE subelusters respeetively. 
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Fig. 3. — Smoothed continuous temperature map from ASCA GIS data of the A1367. The regions used later in the fitting process and 
shown previously in Figure hi are again outlined as is the excluded circular region around NGC 38(52. Pixels with large uncertainties in their 
temperatures ( j,_'^ > 1.5) have been excluded from this map. The SIS data does not extend beyond regions #1 and #7 and the GIS results 
for regions #0 and #8 displayed in Table til and Figure \m are shown for completeness. A scale bar for the temperature is provided for ease of 
reference. 
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ASCA Temperature Profile for A1367 
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Fig. 4. — Fitted temperatures from the two methods with their respective la error bars. Regions #1 through #7 were 5' wide and 16' long 
with the short axis ahgned ta a hne runrung between the SE and NW subclusters. Region #0 and #8 were 5' wide and 21' long with the 
same orientation (see Figure hi or Figure ph. Note the excellent agreement between all of the various results, with the one exception being for 
the GIS data in region #6 with Method A. The large error bar for the SIS data in region #7 is due to the edge of the data. Also GIS results 
for regions #0 and #8 which both lie beyond the SIS frame area are included for completeness. 
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Fig. 5. — Two dimensional temperature map of A1367. The subregions were constructed by dividing the previously defined rectangles into 
three equal parts. Note the numbering scheme was chosen to highlight the association with the previous results. Also included were four 
new regions- #14, #24, #34, and #44- that encompass the cool area to the northeast. The temperatures were determined by simultaneously 
fitting the GIS and SIS data. The peaks of the emission for the SE and NW subclusters are approximately located in subregions #32 and 
#62 respectively. The statistics of all other areas either precluded subdivision-e.g. regions #0 and #8- or simply produced no meaningful 
results even when not subdivided. Again a scale bar for the temperature is provided. Included below the map, for ease of comparison, is a 
plot of the temperatures versus region number with la errorbars. We have drawn lines at 3.2 and 4.2 keV to indicate the weighted averages 
for the two subclusters. 
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Fig. 6. — Fitted abundances from Method A with their respectivej-lcr error bars. Also shown is the abundance used for the results for 
Method B and the continuous map generated from Method A (Figure ph. We note that changing the abundance to 0.4 for Method B had no 
significant effect on the fitted temperatures. 



Table 1 
Observational Data 



Pointing 


a'' 


6'' 


On/Live Time 


b 




Events 








GlS-2/3^ SlS-0 


SlS-1 


aiS-2 


GiS-3 


SlS-0 


SlS-1 


D 


ASCA-SEl 


llM4"52.8^ 


19°41'45.6" 


7,292 6,462 


4,481 


8,482 


8,236 


11,248 


8,053 


o 


ASCA-SE2 


llH4"57.6" 


19°43'19.2" 


9,654 7,452 


6,710 


11,108 


10,881 


13,121 


11,629 


g 


ASCA-NWl 


ll''44™26.4^ 


19°51'57.6" 


10,624 6,944 


5,888 


11,683 


11,664 


10,820 


9,351 


^ 


ASCA-NW2 


ll''44™31.2" 


19°50'24.0" 


8,030 6,490 


4,606 


9,077 


8,996 


10,081 


7,561 


r^ 


ROSAT-PSPC 


llM4'"40.8'' 


19°42'36.0" 


16,357 






203,890 







^J2000 

^Units of seconds. 

^On-Time for the GIS-2 and GIS-3 detectors were the same. 



Table 2 
subcluster data j^ 



Subclustcr 


/3 


Re" 


Tc'^ 


Mtotal'' 


Lo.66Mpc 


Mgas^ 


Teff^ 


ML,,,,^ 




0.5" 1.0" 


0.5" 1.0" 


0.5" 1.0" 


SE 

NW 


0.73 ±0.03 
0.66 ±0.05 


0.42 ±0.02 
0.49 ±0.04 


3.2 

4.2 


0.77 2.3 
0.78 2.5 


0.39 
0.29 


0.061 0.22 
0.049 0.20 


2.2 
2.0 


0.53 1.5 
0.38 1.2 



^Mpc 

'^Masses given are within the radius listed and in units of x IQ^'^Mq 



t3 

0) 



0) 



o 



> 

0) 



'xlO'^^ergs s ^ ^ 



Table 3 
Results of Spectral Fits in the Regions 








1 


2 


3 4 
Temperatures'* 


5 


6 


7 


8 




GIS-A 
GIS-B 

SIS-A 
SIS-B 


3.8 ±0.7 


3.2 ±0.3 
3.9 ±0.6 
3.2 ±0.4 
3.4 ±0.7 


3.2 ±0.2 
3.2 ±0.3 
3.1 ±0.2 
3.4 ±0.5 


3.3 ±0.1 
3.2 ±0.1 
3.1 ±0.2 
3.0 ±0.2 


3.6 ±0.2 

3.7 ±0.2 

4.0+S1 
3.9 ±0.3 


3.8 ±0.3 
3.7 ±0.3 

4.4 ±0.4 


5.0 ±0.5 
4.2 ±0.4 

4.1 ±0.4 
4.1 ±0.5 


/I «+0-8 
4-O-0.6 

4.4 ±0.5 

^ 7+1.6 
^•'-1.0 

4.3 ±0.6 


/I S+l-4 
4-O-0.9 

4.4 ±0.6 


d 

o 
p 

1 










Abundances'^ 










E- 


GIS-A 
SIS-A 


n 0/1 +0-30 


0.23+J>fr 


i^i+oi" 

"•^^-0.09 

o.3ot°:l? 


^•^'-0.07 

S2+°°^ 


^•-■^^-0.08 

14+0-07 

'^•-'^^-0.07 


0.09+ro^ 


0.191°:!^ 
0.08+°;JI 


n AA+0.26 
U-44_o.22 
17+0-18 


n O7+0.36 
"•^'-0.27 





^Temperatures are in keV. The errors quoted are la. 

'^Abundances are relative to solar -with j^''° _ = 4.68 x 10^^ assumed. The errors quoted are la. 



^Total 



